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ABSTRACT 

The frequently observed association between giant radio halos and merging galaxy clusters has 
driven present theoretical models of non- thermal emission from galaxy clusters, which are based on 
the idea that the energy dissipated during cluster-cluster mergers could power the formation of radio 
halos. To quantitatively test the merger-halo connection we present the first statistical study based 
on deep radio data and X-ray observations of a complete X-ray selected sample of galaxy clusters with 
X-ray luminosity > 5 x 10*'' erg/s and redshift 0.2 < z < 0.32. Using several methods to characterize 
cluster substructures, namely the power ratios, centroid shift and X-ray brightness concentration 
parameter, we show that clusters with and without radio halo can be quantitatively differentiated 
in terms of their dynamical properties. In particular, we confirm that radio halos are associated to 
dynamically disturbed clusters and cluster without radio halo are more "relaxed" , with only a couple 
of exceptions where a disturbed cluster does not exhibit a halo. 

Subject headings: galaxies: clusters: general — radiation mechanisms: non-thermal — radio contin- 
uum: general — X-rays: galaxies: clusters 



1. INTRODUCTION 

Radio and X-ray observations of galaxy clusters prove 
that thermal and non-thermal components coexist in the 
intracluster medium (ICM). While X-ray observations re- 
veal thermal emission from diffuse hot gas, radio obser- 
vations of an increasing number of massive galaxy clus- 
ters unveil the presence of ultra-relativistic particles and 
magnetic fields through the detection of diffuse, giant 
Mpc-scalc synchrotron radio halos (RHs) and radio relics 
{e.g., Ferrari et al. 2008; Cassano 2009 for review). RHs 
are the most spectacular evidence of non-thermal com- 
ponents in the ICM. They are giant radio sources located 
in the cluster central regions, with spatial extent similar 
to that of the hot ICM. 

There is collective evidence in the literature that RHs 
are found in clusters with significant substructure in 
the X-ray images, as well as complex gas temperature 
distribution, which are signatures of cluster mergers 
{e.g., Schuecker et al. 2001; Govoni et a. 2004; Marke- 
vitch & Vikhhnin 2001; Venturi et al. 2008). In par- 
ticular, in a seminal paper Buote (2001) provided the 
first quantitative comparison of the dynamical states of 
clusters with RH discovering a correlation between the 
RH luminosity at 1.4 GHz and the magnitude of the 
dipole power ratio Pi/Pq. The RH- merger connection 
suggests that the gravitational process of cluster forma- 
tion may provide the energy to generate the non-thermal 
components in clusters through the acceleration of high- 
energy particles via shocks and turbulence {e.g., Sarazin 
2004; Brunetti et al. 2009). The discovery of RHs with 
very steep spectrum supports the scenario of particle re- 
acceleration by merger-driven turbulence {e.g., Brunetti 
et al. 2008). 

Recently, deep radio observations of a complete sam- 
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pie of galaxy clusters have been carried out as part of 
the Giant Metrewave Radio Telescope ( GMRT) RH Sur- 
vey (Venturi et al. 2007; 2008). These observations 
confirmed that diffuse cluster-scale radio emission is not 
ubiquitous in clusters: only 30% of the X-ray luminous 
{Lx{0.l - 2.4 keV) > 5 x 10"*" erg/s) clusters host a 
RH. Most importantly, these observations allow to sep- 
arate RH clusters from clusters without RH, showing a 
bimodal distribution of these clusters in the 1.4 GHz ra- 
dio power (Pi. 4) versus X-ray luminosity {Lx) diagram 
(Brunetti et al. 2007): RHs trace the well known cor- 
relation between Pi. 4 and Lx, while the upper limits to 
the radio luminosity of clusters with no-RH lie about one 
order of magnitude below that correlation (see e.g., Fig.4 
in Brunetti et al. 2007). The reason for this separation 
is expected to lie in the mechanism responsible for the 
origin of radio emitting electrons (Brunetti et al. 2009 
and references therein). In this Letter, we will show for 
the first time that clusters with RH and clusters without 
RH can be quantitatively differentiated also according 
to their dynamical status. We will use archival Chandra 
data of clusters in the GMRT RH Survey and charac- 
terize cluster substructure in the X-ray images adopting 
different methods. 

A ACDM cosmology {Ho = 70kms-^Mpc-\ = 
0.3, Oa = 0.7) is adopted. 

2. THE SAMPLE AND DATA PREPARATION 

The GMRT RH Survey (Venturi et al. 2007, 2008) is a 
deep, pointed radio survey of clusters selected from the 
ROSAT -ESO Flux Limited X ray (REFLEX; Boringher 
et al. 2004) and extended ROSAT Brightest Cluster 
Sample (eBCS; Ebehiig et al. 1998, 2000) catalogs. 
These two catalogs have almost the same flux limit in the 
0.1 - 2.4 keV band ( > 3 • lO^'^ergs"^ cm'^) and their 
combination yields an homogeneous, flux-limited sample 
of clusters. The GMi?r sample consists of 50 galaxy clus- 
ters with z=0.2 — 0.4, X-ray luminosity Lx > 5 x 10** 
erg/s and declination —30° < S < 60°. 34 clusters in 
the sample had no high sensitivity radio information and 
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were observed with the GMRT at 610 MHz. With the 
above selection criteria the sample is X-ray luminosity 
limited up to 2; ~ 0.25 and X-ray flux hmited at higher 
redshifl|3 see Fig. 1 and 2 in Cassano et al. 2008). Re- 
cently, we have undertaken an extension of the GMRT 
RH survey by considering all clusters in the REFLEX 
and eBCs catalogs with 5 > 60° and with the same z 
and Lx selection. This extension leads to a sample of 67 
galaxy clusters which we refer to as the extended GMRT 
cluster sample. While the radio campaign is ongoing, 
data for 3 clusters of the extended sample are available 
and are considered in the present paper. 

For all clusters in the extended GMRT sample with the 
radio data at hand, we searched in the Ghandra archive 
and found information for a sub-sample of 35 galaxy clus- 
ters. We also required the clusters to have at least 2000 
ACIS-S or ACIS-I counts in the 0.5-2 keV band inside 
an aperture of 500 kpc (see below) in order to produce 
images sufhcient to study the cluster morphological prop- 
erties. Three of the 35 clusters do not match the require- 
ment and thus are not included in our analysis. Clusters 
without radio or X-ray data are unlikely to be selected 
by having/not having a RH or being/not being a merger, 
so the fact that they are omitted should not affect our 
conclusions. Our final sample consists of 32 galaxy clus- 
ters with z=0.2 — 0.4, Lx > 5 x 10''^ erg/s all with radio 
[GMRT and Very Large Array, VLA) and X-ray {Chan- 
dra) data. Only three clusters have z > 0.32. Indeed 
it has been shown that for z > 0.32 the GMRT sam- 
ple is incomplete (Cassano et al. 2008). We will thus 
present results obtained by including and excluding the 
three clusters at z > 0.32. 

All X-ray images have been produced in a standard 
manner using CI AO 4.1.2 (with calibration files from the 
CALDB 4.1.3) in the 0.5-2 keV band. They have been 
renormalized by the exposure maps to maintain unit of 
counts per pixel and preserve a proper application of the 
routines for the estimates of the X-ray morphology pa- 
rameters. We visually inspected each image to remove 
point sources and any residual with low exposure around 
the CCD gaps. 

In this Letter we study the cluster substructure on the 
RH scale analyzing the surface brightness inside an aper- 
ture radius of 500 kpc, since we are interested in the 
cluster dynamical properties on the scales where the en- 
ergy is most likely dissipated. Indeed, studies showed a 
point-to-point correlation between the radio and X-ray 
brightness {e.g., Govoni et al. 2001) and between the 
RH properties and cluster properties (mass, velocity dis- 
persion) calculated within the halo region (Cassano et 
al. 2007). The choice of 500 kpc provides a first natu- 
ral approach, and has also the advantage that it allows 
us to sample both lower and higher redshift clusters in 
the sample with adequate sensitivity. Furthermore, we 
expect a small variation (of about ^ 1.5) in i?A (the 
radius defined as that enclosing a region with an over- 
density A — 200, 500, etc. with respect to the critical 
density at the cluster redshift) among clusters in our sam- 
ple, because they are characterized by very similar X-ray 
luminosity and redshift. This implies that our results 
should not significantly change by considering a radius 



This implies a minimum Lx 
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that takes into account the variation of the cluster ther- 
mal properties. We will explore the morphological esti- 
mators based on different aperture size in a forthcoming 
paper (Cassano et al., in prep.). 

3. MORPHOLOGICAL ESTIMATORS 

The superb angular resolution of Ghandra allows to dis- 
criminate between mergers and relaxed clusters even by 
simple visual inspection. To provide a more quantitative 
measure of the degree of the cluster disturbance, we use 
three methods: power ratios {e.g., Buote & Tsai 1995; 
Jeltema et al. 2005; Ventimiglia et al. 2008; Bohringer 
et al. 2010); the emission centroid shift {e.g., Mohr et 
al. 1993; Poole et al. 2006, O'Hara et al. 2006; Ven- 
timiglia et al. 2008, Maughan et al. 2008, Bohringer et 
al. 2010) and the surface brightness concentration pa- 
rameter {e.g., Santos et al 2008). In the following, we 
briefly describe these methods. 

i) The power ratio method is motivated by the idea that 
the X-ray surface brightness could represent the pro- 
jected mass distribution of the cluster. The power ra- 
tio is a multipole decomposition of the two-dimensional, 
projected mass distribution inside a given aperture Rap- 
Power ratios are usually defined as: 



Po = [ao ln{Rap)] 



(1) 



where aq — S{< Rap) is the total intensity inside the 
aperture radius, and 
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where the moments am and bm are given by: 



and 
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where S{x) is the X-ray surface brightness. Here we will 
make use only of the P^/Pq, that is the lowest power 
ratio moment providing a clear substructure measure 
{e.g., Bohringer et al. 2010). 

ii) Poole et al. (2006) used numerical simulation of clus- 
ter mergers and found that the centroid shift method was 
very sensitive to the dynamical state of the cluster. Fol- 
lowing the method of Poole et al. (2006) and Maughan 
et al. (2008), the centroid shift is computed in a series 
of circular apertures centered on the cluster X-ray peak. 
The radius of the apertures was decreased in steps of 5% 
from Rap — 500 kpc to O.ObRap, and the centroid shift, 
w, was defined as the standard deviation of the projected 
separation between the peak and the centroid in unit of 

Rapi as . 
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where A^ is the distance between the X-ray peak and the 
centroid of the ith aperture. 
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Fig. 1. — a): concentration parameter c vs. centroid shift w 
(see text for details); b): w vs. the power ratio P3/P0', c) P3/P0 
vs. c. Symbols are: RH (red filled dots), no-RH (black open 
dots), mini-halos (blue open dots), z > 0.32 (gray filled dots). 
Vertical and horizontal dashed lines are: c = 0.2, w = 0.012 and 
P3/P0 = 1.2 X 10-7. 

Hi) Following Santos et al. (2008), we made use of the 



concentration parameter, c, defined as the ratio of the 
peak over the ambient surface brightness, S", as: 



S{r < 100 kpc) 
S{< 500 kpc) 



(6) 



The concentration parameter has been used in literature 
for a first identification of cool core clusters in those cases 
where a spatially resolved spectroscopic analysis was not 
possible {e.g., in the case of high redshift clusters; Santos 
et al. 2008). We use the concentration parameter to dif- 
ferentiate galaxy clusters with a compact core (i.e., core 
not disrupted from a recent merger event) from cluster 
with a spread distribution of gas in the core (i.e., core 
disturbed from a recent merger episode). 

It is important to note that among the presented meth- 
ods, the power ratio and the centroid shift methods are 
less sensitive to the presence of substructures (and thus 
mergers) along the line of sight, while the concentration 
parameter is in principle not affected by these projection 
effects. 

4. RESULTS 

We show the results of the substructure analysis in 
Figjl] In Figjl^) we report the distribution of the 32 clus- 
ters m the (c, w) plane. There is a clear anti-correlation 
between the two parameters. Most importantly, RH clus- 
ters (red filled dots) can be well separated from clusters 
without RH (black open dots) and clusters with mini- 
halofP] (blue open dots). In particular, as a reference 
case, II we consider the median value of each parameter, 
w = 0.012 and c — 0.2 (horizontal and vertical lines), 
the sample can be well separated between RH and no- 
RH clusters: no cluster with RH is found in the upper 
region {w < 0.012 and c > 0.2), while the fraction of 
clusters with RH in the lower region {w > 0.012 and 
c < 0.2) is 73-78 % (including or excluding the cluster 
at z > 0.32, respectively). This confirms the hypothesis 
that RH arc located in dynamically disturbed systems. 
We note also that clusters with mini-halos lie in the up- 
per region {w < 0.012 and c > 0.2), supporting the con- 
nection between radio mini-halos and cluster cool cores 
(e.g., Gitti et al. 2002). 

In Fig[T]3) we report the distribution of the 32 clusters 
in the (i(;,P3/Po) plane. We find a clear correlation be- 
tween the two parameters. A similar trend was recently 
found also by Bohringer et al. (2010). Most importantly, 
we find that all RH clusters (the color code is the same 
as above) are located in the region of higher values of the 
parameters w and P^/Pq. The position of clusters with 
mini-halos is consistent with being in the more relaxed 
systems. In this plane the horizontal and vertical lines 
(w — 0.012 and P3/P0 = 1-2 x 10^*") are also the medians 
of each parameter. 

For completeness, in Fig.jlj:) we report the distribution 
of clusters in the (c, P3/P0) plane. There is again a clear 
separation between RH and no-RH clusters, with RH 
clusters located in the region of dynamically disturbed 
systems, with higher values of P3/P0 and lower values 
of c. We note that c is the parameter that provides the 

^ Radio mini-halos are diffuse synchrotron emission on smaller 
scales (e.g., 200-500 kpc) extending around powerful radio galaxies 
at the center of some cool core clusters {e.g., Ferrari et al. 2008). 
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best separation between RH and no-RH clusters, indeed 
no RH is found in clusters with c > 0.2. 

All diagrams provide strong evidence that RHs form 
in dynamically disturbed clusters, while clusters with no 
evidence of diffuse synchrotron emission on Mpc scales 
are more relaxed systems. To test quantitatively this re- 
sult, we run Monte Carlo simulations in the (w, c) plane. 
We randomly distribute the 12 RH clusters among the 
32 clusters of the sample, and count the number of RHs 
falling in the uppcr-lcft quadrant of Fig[I]i) (those with 
w < 0.012 and c > 0.2). We repeat the procedure 10^ 
times and find that only in 3-4 cases no RH is found in the 
upper- left quadrant; this allows us to conclude that the 
observed distribution differs from a random {i.e., inde- 
pendent of cluster dynamics) distribution at more than 
4(7. This proves that our result is statistically signif- 
icant and shows, for the first time, that the separation 
between RH and no-RH clusters has a corresponding sep- 
aration in terms of dynamical properties of host clusters. 
We note that there are 4 outliers in Figjl] (3 if we do 
not consider the cluster at z > 0.32), i.e., clusters that 
are dynamically disturbed but that do not host a RH. 
These clusters deserve further investigation. However, 
their presence in the region of RHs is not surprising in 
the framework of models that explain the cluster-scale 
synchrotron emission with merger-driven turbulence and 
shocks. These models predict that RH should be main- 
tained for a typical lifetime of ~ 1 Gyr (see Brunetti et 
al. 2009) which is of the same order as the merger time- 
scale (during which the cluster would appear disturbed) 
implying that a large fraction of massive and merging 
clusters should host a RH. Most importantly, turbulent 
re-acceleration models predict a cut-off in the spectra 
of RHs at the frequency Vc that is determined by the 
fraction of turbulent energy converted into electron re- 
acceleration. The cut-off makes the observations of RHs 
difficult ai V > h'c- In disturbed clusters with relatively 
smaller masses {My < 2 x 10^^ Mq) and at higher red- 
shifts {z > 0.4 — 0.5) the cutoff frequency can be lower 
(Cassano et al. 2010, Cassano 2010). In line with this 
scenario, 3 out of the 4 outliers have X-ray luminosity 
at the lower boundary of our selection, Lx < 8 x 10^** 

erg/sec {M < 1.9 x 10^^ Mq), and the other is the high- 
est redshift cluster of the sample {z ~ 0.42). They may 



still have RHs, but need lower frequency observations to 
detect them. 

5. CONCLUSIONS 

We used a statistical sample of 32 galaxy clusters with 
Lx{0.1 - 2.4keV) > 5 x 10^^ erg/s with radio {GMRT 
and/or VLA) and X-ray {Chandra) observations, to test 
the merger-RH paradigm in galaxy clusters by relating 
the dynamical state as seen from X-rays to the presence 
of a halo. We adopted three main methods of X-ray 
substructure characterization: the power ratio, P3/P0 
{e.g., Buote & Tsai 1995; Bohringer et al. 2010), the 
centroid shift, w {e.g., Mohr et al. 1993) and the X-ray 
brightness concentration parameter, c {e.g., Santos et al. 
2008). We studied the distributions of clusters in the 
w — c, P^/Pq — w and P^/Pq — c diagrams (see FigfTJ . As 
expected, we found anti correlation between w anoc and 
between c and P3/P0 (clusters with the most compact 
cores are also less disturbed) and a correlation between w 
and P3/P0. RH and no-RH clusters are clearly separated 
in all three diagnostic diagrams, with RHs located in 
more disturbed systems. In particular, the median value 
of each parameter {w ~ 0.012, c ~ 0.2 and P3/P0 ~ 1.2 x 
10^^) splits the sample in RH and no-RH clusters. We 
find no RH cluster in the regions selected by w < 0.012, 
c > 0.2 and P3/P0 < 1.2 x lO""^ (also, no RH is found in 
the region constrained just by c > 0.2), while the fraction 
of RHs increases to ^ 73 — 78% in the regions selected by 
w > 0.012, c < 0.2 and P3/P0 > 1-2 x lO"''. By means of 
Monte Carlo simulations we showed that the probability 
to get such segregation between RH and no-RH clusters 
by chance is of the order of 3 — 4 x 10^^ . This established 
for the first time in a statistical manner the connection 
between RH and cluster mergers. We note also that radio 
mini-halos are located in relaxed clusters (characterized 
by high values of c) supporting the connection between 
radio mini-halos and cluster cool cores. 
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